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CLAIMS 

1. Method of processing digital images comprising 
one or more objects to be quantified, the said method 
comprising the following main stages: 
5 - normalization of the digital images; 

quantization of the images to one bit, further 
comprising at least one of the following stages: 

calculating, from the said images quantized to 
one bit, the perimeter, area and/or fractal dimension of 
10 the said one or more objects to be quantified; 

reconstructing, from the said images quantized 
to one bit, a 3D-image of the said one or more objects 
to be quantified, and/or 

calculating, from the said normalized images, 
15 the fractal dimension of the overall image. 

2. Method according to claim 1, the said method 
comprising a stage of image's normalization (NORM) which 
comprises the following steps: 

1) dividing the image into quadrants; 
20 2) calculating the mean value of intensity of the 

pixels belonging to each quadrant; 

3) calculating the mean value of intensity for 
all the quadrants as a mean of the calculated means of 
step 2) ; 

25 4) setting for each quadrant the mean value of 
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intensity calculated according to step 3) by performing 
one of adding or subtracting a same intensity value to 
each pixel inside a quadrant in order to maintain the 
original Aintensity among the pixels inside a same 
5 quadrant ; 

5) reiterating steps 1) to 4) up to a preset 
quadrant side length. 

3 . Method according to claim 2 , wherein the said 
preset quadrant side length is approximately half length 

10 of the minor side of the said one or more objects to be 
quantified. 

4. Method according to any one of claims from 1 
to 3, wherein the said digital image has been acquired 
by a confocal microscopy. 

15 5. Method according to claim 4 # wherein the said 

confocal microscopy is a Laser Scanning Confocal 
Microscopy (LSCM) or a Scanning Ophtalmoscopy . 

6. Method according to any one of claims from 1 
to 5, which comprises the following steps: 
20 la) dividing the image into four quadrants; 

2a) calculating the mean value of intensity of the 
pixels belonging to each quadrant; 

3a) calculating the mean value of intensity for 
the four quadrants as a mean of the four calculated 
25 means of step 2a) ; 
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4a) setting for each quadrant the mean value of 
intensity calculated according to step 3a) by performing 
one of adding or subtracting a same intensity value to 
each pixel inside a quadrant in order to maintain the 
5 original A in tensity among the pixels inside a same 
quadrant ; 

5a) determining for each quadrant the max and the 
min values of intensity of the pixels and calculating 
for each pixel an extended intensity value (EI) which 
10 derives from the stretching of the digital values inside 
the range of the possible digital values; 

6a) setting for each pixel the EI p i X ei calculated 
according to step 5a) ; 

7a) reiterating steps la) to 6a) up to a preset 
15 quadrant side length. 

7. Method according to claim 6, wherein the said 
step 5a) of calculating the EI value of the pixels is 
performed by means of the following algorithm: 

Elpixel = (Ipixel ~ Imin) X N / ( Iraax - Imin) 

2 0 wherein Ipixei is the intensity of each pixel of a 

given quadrant, Imin is the min value of intensity of the 
pixel inside the said quadrant, Imax is the max value of 
intensity of the pixel inside the same quadrant and N is 
an integer more than 1 and up to 255, preferably 255. 

25 8. Method according to any one of claims from 1 
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to 5, wherein the said normalization stage comprises: 
lb) dividing the image into quadrants ; 
2b) determining for each quadrant the max and the 
min values of intensity of the pixels and calculating 
5 for each pixel an extended intensity value (EI) which 
derives from the stretching of the digital values inside 
the range of the possible digital values; 

3b) storing the EI P i xe i value for each pixel of each 
quadrant in a data structure; 
10 4b) reiterating steps lb) to 3b) up to a preset 

quadrant side length in order to obtain for each pixel a 
set of intensity values in the data structure; 

5b) calculating for each pixel the mean of the 
intensity values of the set stored in the data structure 
15 and setting the calculated mean value to the respective 
pixel . 

9 . Method according to claim 8 , wherein the said 
step 5a) of calculating the EI value of the pixels is 
performed by means of the following algorithm: 

20 Elpixel = (Ipixel - Imin) X N / ( Imax " Imin) 

wherein I piX ei is the intensity of each pixel of a 
given quadrant, I min is the min value of intensity of the 
pixel inside the said quadrant, Imax is the max value of 
intensity of the pixel inside the same quadrant and N is 
25 an integer more than 1 and up to 255, preferably 255. 
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10. Method according to any one of claims from 1 
to 9, further comprising a stage of image elaboration 
(IMA- EL stage) to quantize the image to "1 bit" . 

11. Method according claim 10, wherein the said 
5 IMA-EL stage comprises the following steps: 

lc) considering a parameter for each pixel; 
2c) comparing said pixel's parameter with a preset 
threshold value or threshold range for said parameter; 

3c) selecting a cluster of active pixels and a 
10 cluster of inactive pixels on the base of said 
comparison, 

wherein said pixel's parameter is preferably 
brightness intensity (black and white images) or digital 
colour value. 

15 12. Method according to any one of claims from 1 

to 11, further comprising a stage of image 
quantification which comprising at least one of the 
following steps: 

- calculating the area A of the object under 
20 examination by counting the number of pixels belonging 

to the cluster of active pixels selected according to 
the previous IMA- EL stage; 

- calculating the perimeter P of the object under 
examination by i) selecting the object contour's pixels, 

25 and ii) applying to such selected pixels a perimeter 
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calculation's algorithm, wherein to each active pixel 
belonging to the object is given a "perimeter value", 
which is a function of the position of the active pixels 
adjacent to the pixel under examination, the sum of said 
5 tt perimeter values" being the overall perimeter P of the 
object . 

13 . Method according to any one of claims from 1 
to 12, further comprising a stage of object's sorting 
(SORT) for identifying objects made up from 4 -connected 
10 pixels, which includes the following steps: 

Id) scanning of the image quantized to w l bit" 
along a predefined direction on a x, y axis system; 

2d) selecting a first active pixel along said 
direction of scanning, said active pixel being 
15 identified by a first set of x, y values, said first 
active pixel belonging to a first object's image; 

3d) performing on said first selected active pixel 
a search routine in the positions next to said selected 
pixel on the direction's line; 
20 4d) iterating step 3d) until an inactive pixel is 

found ; 

5d) assigning to each active pixel selected 
according to such steps 3d) and 4d) a set of x, y 
values, saving them in the storing means of the 
25 processing system (7) and switching said pixels from 
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active to inactive in the object's image; 

6d) evaluating for each pixel selected according 
to steps 3d) , 4d) and 5d) the two next pixels in the 
direction orthogonal to the said scanning direction and 
5 selecting the active pixels; 

7d) performing, for each of said active pixels 
selected according to step 6d) , the routine of steps 3d) 
to 5d) ; 

8d) iterating steps 6d) and 7d) until all of the 
10 connected pixels belonging to the same object have been 
saved; 

9d) repeating steps Id) and 2d) until a first 
active pixel of a further object's image is found; 

lOd) repeating steps 3d) to 9d) until the whole 
15 image has been scanned. 

14. Method according to claim 13, wherein said 
predefined direction in step Id) is from left to right 
starting from top to bottom. 

15. Method according to claim 13 or claim 14, 
20 wherein the stage of object's sorting according to steps 

Id) to lOd) is performed for also identifying objects 
made up from 8-connected pixels, in said stage the step 
6d) being modified as follows: 

6d) evaluating for each pixel selected according 
25 to steps 3d) , 4d) and 5d) the two next pixels in the 
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direction orthogonal to the said scanning direction and 
the two pixels adjacent to each of these latter pixels 
on the parallel line adjacent to the direction's line 
and selecting the active pixels. 
5 16. Method according to any one of claims from 13 

to 15, further comprising at least one of the following 
steps : 

le) calculating the area of each object identified 
according to the SORT stage by counting the number of 
10 pixels belonging to said object's image and multiplying 
it for the area of each pixel; and/or 

2e) counting the number of objects and calculating 
its density; and/or 

3e) calculating the mean area of the objects by 
15 adding the areas calculated according to step le) of all 
the objects sorted and dividing the total area by the 
number of objects obtained according to step 2e) . 

17 . Method according to any one of claims from 1 
to 16, further comprising a step of calculating a 
2 0 parameter (w) indicating the degree of "rugosity" of the 
selected object, the said (w) parameter being preferably 
calculated by means of the following algorithm: 

w = Pf/ R 

/l^jAf.n 

wherein Pf is the perimeter, Af is the area of the 
25 object and R is the "roundness coefficient" of the 



WO 2005/006255 



PCT/IB2003/002703 



41 

object; wherein R is on its turn calculated with the 

following algorithm 

R = Pe/ 

/ l^Aejt 

wherein Pe is the perimeter of the ellipse in which 
5 the measured object is inscribed and Ae its area* 

18, Method according to any one of claims from 1 
to 17 , further comprising a stage of dimensional 
calculation (DIM-CLC) for calculating the fractal 
dimensions of perimeter and area of the observed 
10 objects, wherein said fractal dimension of the perimeter 
(Dp) and said fractal dimension of the area <D A ) are 
determined according to the following steps: 

a) dividing the image of the object into a 
plurality of grids of boxes having a side length e, in 

15 which s varies from a first value substantially 
corresponding to the side of the box in which said 
object is inscribed and a predefined value which is a 
fraction of said first value, 

b) calculating a value of a logarithmic 
20 function of N(e) , in which N(e) is the number of boxes 

necessary to completely cover the perimeter (P) or the 
area (A) , respectively, of the object and of a 
logarithmic function of 1/s for each e value of step a) , 
thus obtaining a first set of values for said 
25 logarithmic function of N(s) and a second set of values 



WO 2005/006255 



PCT/IB2003/002703 



42 

for said logarithmic function of l/e, 

c) calculating the fractal dimensions (D P ) or 

(D A ) as the slope of the straight line interpolating 
said first set of values for said logarithmic function 
5 of N(e) for the perimeter (P) or the area (A) , 
respectively, versus said second set of values of step 
b) . 

19. Method according to any one of claims from 1 
to 18, further comprising a stage of surface 
10 quantification (S -QUANT) performed on the image 
normalized according to the NORM stage, the said stage 
comprising the following steps: 

If) dividing the image in a x, y bidimensional 
.mesh with n x n boxes of side 2; 
15 2 f ) dividing the 0-256 grey scale into n 

subregions having each a 256/n value; 

3f) calculating for each box of the x, y 
bidimensioanl mesh the min and max value of the pixels 
contained therein and of the pixels that contour the 
20 box; 

4f ) calculating how many subregions of 256/n value 
are included between the min and max values of the 
pixels of each box; 

5f) calculating the number N(l) of tridimensional 
25 boxes of side 1 that intercepts the image's surface as a 
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sum of the subregions of all the boxes calculated 
according to step 4f ) ; 

6f) reiterating steps If) to 5f) with a side length 
1' less than 1; 

5 7f ) by repeating step 6f ) , generating a first set 

of values of a logarithmic function of l/l and a second 
set of values of a logarithmic function of N(l) ; 

8f) calculating the fractal dimension of the 
image's surface as the slope of the straight line 
10 interpolating said first set of values versus said 
second set of values of step 7f ) . 

20. Method according to any one of claims from 1 
to 19, further comprising a stage of 3D- reconstruction 
(3D-R) performed on the image subjected to the IMA- EL 
15 stage, the said 3D-R stage comprising the following 
steps : 

Ig) overlapping each image with the subsequent 
image along the z axis; 

2g) minimizing the difference of brightness and/or 
2 0 colour intensity between overlapping pixels by shifting 
along the x axis and/or the y axis an image with respect 
to each other; 

3g) repeating steps lg) and 2g) for each pair of 
ad j acent images . 
25 21. Method according to claim 20, further 
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comprising a stage of object counting ( O- COUNT) , which 
comprises the following steps: 

lh) scanning of the 3D-image quantized to *l bit" 
along a predefined direction on a x, y axis system; 
5 2h) selecting a first active pixel along said 

direction of scanning, said active pixel being 
identified by a first set of x, y values, said first 
active pixel belonging to a first object's image; 

3h) performing on said* first selected active pixel 
10 a search routine in the positions next to said selected 
pixel on the direction's line; 

4h) iterating step 3h) until an inactive pixel is 
found; 

5h) assigning to each active pixel selected 
15 according to such steps 3h) and 4h) a set of x, y 
values, saving them in the storing means of the 
processing system 7 (all of such pixels will have the 
same y value and x values in progressive order) and 
switching said pixels from active to inactive in the 
20 object's image; 

6h) evaluating for each pixel selected according to 
steps 3h) , 4h) and 5h) the two next pixels in the 
coplanar direction orthogonal to the said scanning 
direction and the two next pixels along the z axis, in 
25 the directions +z and -z, and selecting the active 
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7h) performing, for each of said active pixels 
selected according to step 6h) , the routine of steps 3h) 
to 5h) ; 

5 8h) iterating steps 6h) and 7h) until all of the 

connected pixels belonging to the same object have been 
saved ; 

9h) repeating steps lh) and 2h) until a first 
active pixel of a further object's image is found; 
10 lOh) repeating steps 3h) to 9h) until the whole 

image has been scanned; 

llh) counting of the number of the objects sorted 
according to steps lh) to lOh) . 

22. Method according to claim 21, wherein the said 
15 predefined direction in step lh) is from left to right 

starting from top to bottom. 

23. Method according to claim 21 or claim 22, for 
sorting also 8 -connected pixel objects, wherein step 6h) 
of the procedure depicted in claim 21 is modified as 

20 follows: 

6h) evaluating for each pixel selected according 
to steps 3h) , 4h) and 5h) the two next pixels in the 
coplanar direction orthogonal to the said scanning 
direction and the two next pixels along the z axis, in 
25 the directions +z and -z, and the two pixels adjacent to 
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each of these pixels on the parallel line adjacent to 
the direction's line and selecting the active pixels. 

24 . Method according to any one of claims from 1 
to 23, further comprising a stage of volume calculation 
5 (V-CLC) which comprises the following steps: 

li) calculating the area of each object in a first 
2D-image corresponding to a first object's section; 

2i) multiplying the area calculated according to 
step li) for the distance between the said first 
10 section's image and the subsequent section's image, 
taken in the z direction of scanning, wherein an image 
of the same object is contained ; 

3i) reiterating steps li) and 2i) for each 
section's image in the order. 
15 25. Method according to claim 24, wherein the 

overall volume of the objects in the examined tissue is 
determined as the sum of the single volumes. 

26. Method according to claim 24 or claim 25, 
wherein the volume is calculated as : 

20 v = \IM(A + a + 4Axi) 

wherein d is the distance between the two sections, 
A is the area of the first object's section and a is the 
area of the second object's section. 

27. A system (1) for acquiring and processing an 
25 image including a confocal scanning microscope (2) , 
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electronic image acquisition means (6) operatively 
connected to said microscope (2) , a processing system 
(7) operatively connected with said confocal scaning 
microscope (2) and said image acquisition means (6) , 
5 said processing system (7) comprising a processing unit 
(CPU) , storing means which include a RAM working memory 
and a hard disk, said processing system (7) running a 
program (PRG) to perform a method according to any one 
of claims from 1 to 26. 

10 28. A software program (PRG) to perform the method 

according to any one of claims from 1 to 26. 

29 . A computer readable support comprising a 
program (PRG) to perform the method according to any one 
of claims from 1 to 26. 

15 30. Use of a system (1) according to claim 27 or 

claim 28 , for performing a method as depicted in any one 
of claims from 1 to 26. 



